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Abstract 



o 

OA ' In this work, we consider homogeneous and isotropic FRW model of the universe, filled with 

interacting dark matter and dark energy. The dark matter is chosen as usual in the form of dust 
while dark energy is holographic in nature with IR cut off at the Ricci's length and it is in the 
form of a perfect fluid with variable equation of state. We have chosen the interaction term of the 
' following two types: (i) a linear combination of the matter density of the two fluids, (ii) a product 

(N ■ of the two matter densities. For both the choices the evolution equations are transformed to an 

autonomous system and the corresponding critical points are analyzed. Finally, for the first choice 
O , of the interaction term the evolution of the ratio of the energy densities has been studied from the 

point of view of the present coincidence problem. 
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^ ! 1 Introduction 

In the last decade there are wide variety of modern cosmological observations namely precision measurements 
of anisotropics in the cosmological microwave background radiation 1 , baryon acoustic oscillation [31 [5] and 
type la supernova [4l[5l[6]. These observations indicate that at present our universe is composed of nearly 25% 
cold dark matter (DM), 70% nonbaryonic unknown matter known as dark energy (DE)[4l O [71 HI [9] and 5% 
of radiation and baryonic matter which ia well understood by the standard models of particles. The natural 
and leading choice of the unknown DE is the cosmological constant(ACZ3M model) which represents a vacuum 
^ ' energy density having constant equation of state u — —1. However, its observed value is far below than the 

esteemation from quantum field theory (known as cosmological constant problem). Also there is no expectation 
why the constant vacuum energy and matter energy densities are precisely of the same order today(coincidence 
problem). Due to these observational [T0| and theoretical [HI [E] probes for the cosmological constant there 
are alternative models for DE (varies from time) in the literature. Scalar field models [13, 14] E] (known as 
quintesence) have attracted special attention compared to the other alternatives [12^. 

The fact that at present DE and DM are dominant sources of the content of the universe, there has been a 
lot of interest in studying coupling in the dark sector components [T6l[T7l[T8lP^[Ml[2T l [22 l [23 l [24 l [27 l [25l[26l 
[28l [29l l30l [3TI l32] . It is partly motivated by the fact that one can only extract information of these components 
through gravitational interaction. Also consideration of interaction is natural in the framework of field theory 
[m [19] . Recently, it has been shown that an appropriate choice of the interaction between DE and DM can 
alleviate the coincidence problem [ ^ [ ^ [ ^ [ ^ 1 ^ . 

To have some inside about the unknown and mysterious nature of DE, many people have suggested that DE 
should be compatible with Holographic principle , namely "the number of relevant degrees of freedom of a system 
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dominated by gravity must vary along with the area of the surface bounding the system" |25) . Such a DE model 
is known as Holographic DE(HDE) model. Further the energy density of any given region should be bound 
by that ascribed to a Schwarzschild black hole(BH) that fills the same volume [2^. Mathematically, we write 

PD < MpL~^, where pa is the DE density, L is the size of the region(or infrared cut off) and Mp — (SttG) ^ 
is the reduced Planck mass. Usually, the DE density is written as 

PD = -^2- (1) 

Here the dimensionless parameter 'c^' takes care of the uncertainties of the theory and for mathematical con- 
venience the factor 3 has been introduced. Due to lack of clear idea there are many choices for the infrared 
cut off of which the most relevant one are the Hubble radius, i.e., L = if~^[3TJ [37] and the Ricci's length, i.e., 



L = i^H + 2H^j [Ml EH IMl ED- The argument behind the choice of Ricci's length as the IR cut-off is that 
it corresponds to the size of the maximal perturbation, leading to the formation of a black hole[32]. Another 
commonly used IR cut-off length is the radius of future event horizon, but it suffers from a severe circularity 
problem. 

In the present paper we consider a cosmological model of Holographic DE (HDE) in the form of a perfect 
fluid interacting with DM in the form of dust. The choice of the IR cut-off is chosen as the Ricci's length, i.e.,|33] 



L=[H + 2H^j ' (2) 

The evolution equations of the model are formulated into an autonomous system and critical points are analysed. 
Explicit solutions are obtained and are analysed asymptotically. 



2 Basic equations for interacting HDE at the Ricci scale 

In the present work, the homogeneous and isotropic FRW universe is assumed to fill up with interacting two 
fluid system-one component is in the form of dust (having energy density p^) known as dark matter (DM) while 
perfect fluid having barotropic equation of state pn = uj£,po, ojd, a variable is the DE component. 
Assuming spatialy flat model, the friedmann equations are (choosing SttG = 1 = c) 

= + (3) 

and 

2H=-Pm-{1+UJD)PD (4) 

and we have the conservation equations 

Pm + 3i?Pm = Q (5) 

PD + iH{l+UJD)pD = -Q (6) 

Here the interaction term Q > Q indicates transfer of energy from DE component to DM sector while opposite 
is the situation for Q < 0. As Q < would worsen the coincidence problem so we choose Q > throughout 
the work. Also validity of the second law of thermodynamics and Le chatelier's principle [TBI IM] support this 
choice of positive Q. It should be noted that we have not included baryonic matter in the interaction due to 
the constraints imposed by local gravity measurements [13l [34] [35] . In the next two sections we shall deal with 
two different choices of interaction term separately, namely, (i) Q = Sb^Hp [p — pm + Pd, the total energy 
density) and (ii) Q = jpmPD {j > 0). The first choice is the special case of the usual one used in the literature 
as a linear combination of the energy densities. The second choice is physically more viable in the sense that 
interaction rate vanishes if one of the densities is zero and increases with each of the densities. Also this choice 
of interaction for HDE models gives the best fit to observations [El [34]. Also it should be noted that the 
constant 7 has the dimension of [L'^AfT] . 

Using the field equations ([3]) and ((H) we have the form ([T]) the expression for the energy density of HDE as 

= y {Pm + (1 - 3a;£)) P_d} (7) 
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Hence the equation of state paranietercan be expressed in terms of the density parameter as 



1 



ud = 



3c2 3n 



D 



Also the deceleration parameter takes the simple form 

a 



= 1 ^ = 1- 



1 



c2 (1 + u) 



(8) 



(9) 



which shows a smooth transition from deceleration to acceleration as universe evolves from the early matter 
dominated era to the late time DE dominated and here u = — = tt 1- 

Pd iir> 



3 Explicit calculations for choice of interaction term 

Case (I) : g = Sb^Hp 

The energy conservation equations can be written explicitly as 



Pm = \/3 {prn + Pd) [h^ PD - (l - b^) Pn 



Pd = -SH [b'^ {p,n + Pd) + Pd + ^^d)] = -^3 {Pm + Pd) 



1 



b ip,n + Pd)+ [1 + ^ T^] PD 



3c' 



(10) 



(11) 



As a consequence the evolution of the density parameter Qd and the ratio of the energy densities u has the 
form 

2nD (1 -^d) , ^ , , .J 



n, 



-H 



u = H 



+ (1 - + 36" 



-^ + u(l + u) + 362(l + u)^ 



(12) 



(13) 



From equation (jl2p we see that the DE density gradually decreases with the evolution of the universe until 
it is in the phantom era (1 + bJD < 0). Thus if we assume the DE density to be sufficiently large at the 
early epochs of the evolution then from equation (jlip p'^ — increases till some intermediate stage and then 
gradually decreases with pm — along the line — = j^p- in the (p,„, pd) plane. Thus we have u 0(1) at 
some intermediate stage of evolution in the neighbourhood of the above line and may be a possible resolution 
of coincidence problem |33[ 134] . Further, the present model of the universe shows a DE dominance at the early 
epochs and subsequently universe evolves with DM as the dominant component and then again it has DE domi- 
nated phase at late time as predicted by observation. Hence the above model is suitable for the present universe. 

Now to formulate an autonomous system we rewrite equation ([T2l) using field equation ([3]) as 



p'd 



3H^ {I + 3b^) + PD 3-- 



and the second field equation, i.e., equation (|4]) can be written as 



(14) 



(15) 



where ' = -§-, x = Ina Thus equations (|14|) and ()15|) form a linear homogeneous autonomous system in 



the phase plane [pD, H^j, having critical point at the origin. For the Jacobi matrix A, Tr{A) — — 7- 
det{A) = ^(6c2- 5-362). 



and 



The nature of the critical points is characterized |36j in the Tablel and the geometrical features are pre- 
sented in fig 1(a) and 1(b). 
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Phase Portrait for /r=0.294 and 6c^<5+3fe^ , , 2 h- 5 

— , — , — , — „ Phase Portrait for h =0.294 and c >min( - , — + - ) 

^^^^ H '7 2 6 




^ LC'^/^"- ' ^'^^^^ ^^^^^ 



Fig. 1(a) . Fig. 1(b) 



Fig. 1(a)- 1(b) represent the variation of pD — H. Though in Fig 1(b) the whole region is not a 
physically valid region but for better understanding about the system and the nature of the critical 
point we have drawn graph for the whole region. 
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Table 1: Nature of the Critical Points for Case-I 



Condition 


Nature of the eigen values 


Type of critical point 


(i) 6c^ < 5 + 36^ 

{a) > min (^f , ^ + i) 


(i) real roots of opposite sign 

(ii) both negative real roots 


(i) Saddle 

(a) Stable nodes. 



From equation (jl3p if Uf be a fixed point, i.e., = then the parameter has the expression 

3 {1 + Uff 

Using this value of 6^ the other fixed point of equation (|13l) can be expressed in terms of u/ as 



2 

7^' 



-(1 + ^/) 

rUf + (1 + U/) 

Further from equation (jl3p . differentiating once we obtain 

du' 2 



du 



+ 1 + 2m + 66'' (1 + u) 



(16) 



(17) 



(18) 



Now if we assume the fixed points Up and Uf to be at the far past and at the far future respectively, then 
Uf < uq ~ 0.45 and < — 0.75 [33], where the suffix '0' stands for the present value and the second 
inequality is obtained from equation (jO]) with the fact that we are at present in an accelerating phase. Then 
from 1181) 



du' 
du 



(2 + c^) 



and 



du!_ 
du 



1 



C2(l+M/) 

2 [uf - 1) 



> 



C2 {uf + 1) 



< 



So the fixed point at far past is an unstable one while the fixed point at far future is a stable one. 
Moreover, integrating equation (1131) we obtain 



(Up + Uf) {uf - Up) 



l + k{d)^' 
k{a)f^ - 1 



(19) 



where k = d = ^and^i={l + 36^) {up - Uf). 

Further, equation (jl3p can be expressed in terms of Up and Uf as 

u = (l + 36^) (u — Up) {u — Uf) . 



(20) 



which clearly shows that ii < between the two fixed points Up and Uf. The continuous decrease of u 
between the two fixed points is shown in figure 2(a) where u ~ 1 near a = 0.8. Hence the coincidence problem 
has some partial solution for the present model, it can not predict Uq ~ 0(1) [33) . The explicit expression for 
the density parameter O/j is given by 



k(d)^ 



1 



and hence using (|S]) we have 



2 
3? 



(1 + 1*/) fc(a)^ - {l + Up) 

l{l + Uf)k{h)^' -{l + Up) 



k{a)^ 



(21) 



(22) 
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Plot for u against a Plot for % against a 
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Fig.2(a) . Fig.2(b) 



Fig. 2(a)-2(b) represent the variation of u and lod against a respectively corresponding to the value 
of the parameters Uf = 0.013, = 0.44, uq = 0.4144 and 6^ = 0.290. 



The variation of W£) over the scale factor is shown in figure 2(b) which shows that we are very close to 
AC DM era in the present epoch. At the two extreme limits the limiting values of are as follows : 



2 

3? 



UJo — > LUoo — 



jl+Up) 

3 

(1+Mp) 



(a ^ 0) 



2 

3? ' 3 V- ' 

Note that as u decreases with the evolution so it/ < Up and hence ujjj also decreases as the universe grows 
up. If we choose mq ~ 1 then the present value of ujd (i-e.,WDo) does not depend on the asymptotic values Up 
and Uf, i.e., 



which is compatible with recent observation, i.e., ujjjq should be very close to —1 (as shown in Fig. 2(b)) if 
~ 0.4 and is closed to the estimated lower bound of c^. Integrating field equation ^ using equation ([2T|) 
and (1221) we have 



(1 + u/)/ca^ - (1 + Up) 



L {1+Uf)k-{1+Up) 

Now combining equation (|2ip and the HDE density has the form 



(23) 



Pd 



[Up - uf) {uo + 1) 



{up - Uo) (a) 



+ (t^o - Uf) (a) 



(24) 
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Table 2: Nature of the Critical Points for Case-II 



Condition 


Nature of the eigen values 


Type of critical point 


(i) 5^ < < i 

(ii) (? >\ 

(iii) <k 

(^^) = 
(.) c2 = i 


(i) pair of complex roots with real part positive. 

(ii) real roots of opposite sign. 

(iii) real roots of same sign. 
(iii) equal non-zero real roots. 
(iii) zero roots. 


(i) unstable spiral. 

(ii) Saddel. 

(iii) Unstable node. 

(iii) degenerate unstable node. 
(iii) degenerate point. 



The above expression for pu contains two terms — the first one is dominant at later epochs when a is large 
while the second term is the dominant one at early phases. 

Case (II) : Q = 7p,„pD 

As before the explicit form of the energy conservation equations are 

Pm. = Pm [iPD - iH] (25) 

and 



PD = -PD [iPni + iH(l + UJD)] 

and hence the ratio of the energy densities has the evolution equation 

(1 + ^.) 



ii = 3Hu ['-/H 
The field equation (U) can be written as 



3Hu 



2 

3? 



H = -- 



3c2(f 



(26) 



(27) 



(28) 



Thus equations (j27p and (j28p from an autonomous system in the (u, _ff)-phase plane. The only critical point 

which is of physical interest is (^-^ — 1, 2^~2^ in the (u, H) phase plane. It is to be noted that the other 

critical points correspond to static model of the universe or a degenerate line |36j representing universe filled up 
with DE only. The critical point is characterized in the Table2 where for the linearized matrix A, p = tr{A) = 



1 



(i - c^), g = detA 



c^) and A= - Aq 



v50 



the equilibrium point in each of the above five cases are presented in figures 3(a)-(e) 



c ) . The geometric nature of 



4 Discussion: 

The paper analyzes the HDE model interacting with DM (in the form of dust). Here the IR cut off is chosen 
at the Ricci's length with the justification that it corresponds to the size of the maximal perturbation corre- 
sponding to formation of a black hole. The interaction between the two fiuids is either a linear combination 
or in product form of the two energy densities of which the product form is physically more variable one. In 
both the cases, the evolution equations are transformed to an autonomous system for which the nature of the 
critical points are presented in tabular form and are graphically analyzed. Finally, the coincidence problem is 
discussed by studying the evolution of the energy ratio for the first case only. 



7 



1 1 

Phase Portrait for y=0.232 and — <c < - 
50 2 




/ i i i \ \ \ \ \ \\ \^ 

y n n t I i \ \ \ \ \ 
/////// M \ \ \ - 
, //////// i i \ \ 



Phase Portrait for 7=0.232 and c >- 




\ ^--^ — — '^V- 

\ \ \ \ \ / /, - 

\ \ \\\\\\\ \ \ //^>^ 



Fig.3(a) 



Fig.3(b) 
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Phase Portrait for r=0.232 and c^= 
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